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Stirling engine is an external combustion engine that uses any external heat source to generate 
mechanical power which operates at closed cycles. These engines are good choices for using in power 
generation systems; because these engines present a reasonable theoretical efficiency which can be 
closer to the Carnot efficiency, comparing with other reciprocating thermal engines. Hence, many studies 
have been conducted on Stirling engines and the third order thermodynamic analysis is one of them. In 
this study, multi-objective optimization with four decision variables including the temperature of heat 
source, stroke, mean effective pressure, and the engine frequency were applied in order to increase the 
efficiency and output power and reduce the pressure drop. Three decision-making procedures were 
applied to optimize the answers from the results. At last, the applied methods were compared with 
the results obtained of one experimental work and a good agreement was observed. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

According to the growing enhancement of fossil fuels consump¬ 
tion and limited energy resources in the world, the environmental 
pollution, and global warming, the attentions have been focused on 
the field of energy, especially renewable energy, in recent years. 
Solar energy is one of the most important sources of renewable 
energies, and Stirling engine is one of the functional techniques 
that uses solar energy [1,2]. Stirling engine is an external heat 
engine that works at closed cycles. The heat from the external heat 
source is transferred to working fluid and is converted to work via 
the expansion of the working fluid within the cylinder. Since Stir¬ 
ling engine is an external heat engine, a variety of thermal 
resources like solar energy, nuclear energy, the heat produced from 
the combustion of fossil fuels and waste heat from industrial pro¬ 
cesses can be used as the external heat source [3], High thermal 
efficiency, minimal pollution, reliability, maximum utilization 
and clean combustion are the major features of Stirling engines. 
Due to the potential of Stirling engine which was mentioned above, 
this engine is a good alternative to the conventional power gener¬ 
ation systems. Hence, many studies have been carried out on the 
modeling of Stirling engines of which, some are mentioned in the 
following sections. 


* Corresponding author. Tel.: +98 9121947510; fax: +98 21 88617087. 
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The first analysis on the performance of Stirling engines was 
presented by Schmidt in 1871. The Schmitt’s theory is an effective 
tool in the design of Stirling engines 4]. The assumptions that 
Schmidt considered in his analysis included: There is no tempera¬ 
ture gradient in heat exchangers, the surface temperature of the 
cylinder and the piston is constant, the mass of the fluid is constant 
which implies that there is no leakage, the perfect gas equation of 
state is applied, the volume changes are sinusoidal, and the speed 
of working fluid within the engine is considered constant. One of 
the aforementioned assumptions was that the gas behaves isother¬ 
mal in the expansion and compression spaces. In practice, this 
assumption is not true, especially at high speeds. At the higher 
speeds, the expansion and compression processes are closer to adi¬ 
abatic process. Based on this fact, Finkelstein presented adiabatic 
model (the second order model) for creating a more realistic anal¬ 
ysis of Stirling engine [5], In adiabatic model, the engine is divided 
to five parts (cooler, heater, regenerator, expansion and compres¬ 
sion spaces). Each part of the engine is considered as a control vol¬ 
ume and the conservation laws of mass and energy are analyzed 
for each part [5]. 

Another analysis of Stirling engine is the third order analysis. In 
the third order model, the equations of mass, momentum and 
energy of engine are written in the form of partial differential 
equations for several control volumes. Some examples of such 
models were investigated by Dyson et al. [4] and Anderson [6], 
Firstly, this model was applied by Feurer for designing Stirling 
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Nomenclature 



d 

displacer stroke (cm) 

Vf 

gas velocity (m/s) 

r c 

crank radius (m) 

P 

gas density (kg/m 3 ) 

A f 

gas flow area (m 2 ) 

n 

normal vector on control-volume surface 

f 

engine frequency (Hz) 

X 

control-volume length (m) 

n 

engine efficiency 

Aw 

wall area of control-volume (m 2 ) 

p 

engine output power (W) 

s 

connecting surfaces between adjacent control-volumes 

Or 

total heat rate rejected from the engine (W) 


(m 2 ) 

T c 

heat rejection temperature (K) 

V 

volume of control-volume (m 3 ) 

T e 

heat absorption temperature (K) 

T 

gas temperature inside control-volume (I<) 

Re 

Reynolds number of control-volume gas 

P 

gas pressure inside control-volume (Pa) 

R6 w 

Reynolds number of control-volume gas 

if 

gas shear stress from control-volume wall (Pa) 

p 

dynamic viscosity of control-volume gas (Pa s) 

u 

gas specific internal energy (J/kg) 

Nu 

Nusselt number 

Qc 

net absorbed heat of control-volume gas (J) 

k 

control volume index 

Qe 

net injected heat to control-volume (J) 

1< 

Working gas conductivity (W/m K) 

m w 

mass of control-volume wall (kg) 

t 

time (s) 

c w 

specific heat of control-volume wall (kj/kg) 

d i+ 

distance from the positive ideal point 

T w 

temperature of control-volume wall (K) 

dr 

distance from the negative ideal point 

r 

working gas constant (kj/kg K) 

c v 

specific heat at constant volume (kj/kg K) 

T h 

temperature of the heat source (K) 

L 

length of connecting rod (m) 

Vde 

expansion dead volume 

e 

eccentricity 

Ad 

sectional area of the displacer piston (m 2 ) 

Vole 

compression dead volume 

A p 

sectional area of the power piston (m 2 ) 

D h 

hydraulic diameter (m) 

dP 

Pressure drop (bar) 


engine [7 . The third order analysis is used for the optimization 
methods especially for evaluating the efficiency, output power, 
and the behavior of the working fluid flow inside the engine. So 
far, several algorithms have been used for optimization of Stirling 
engines. Each of these algorithms is attended on the part of the 
predicted results by a mathematical model. By recording the 
changes of the results, the orientation of the process optimization 
and the optimal values of the design variables were determined. 
For example, Carlsen by using two variables power and thermal 
efficiency that was obtained from a mathematical model increased 
the efficiency of Stirling engine up to 23% [8], Chen used genetic 
algorithms to optimize a solar Stirling engine and the results of 
the optimization were implemented on the engine [9], 

Kraitong and Mahkamov carried out theoretical investigations 
on the determination of optimal design parameters of a low tem¬ 
perature difference (LTD) solar Stirling engine using the optimiza¬ 
tion method based on Genetic algorithms [10], For considering all 
the aforementioned issues, we followed the third-order model and 
studied multi-objective functions including the pressure drop 
(bar), the output power (kW) and the efficiency of the considered 
Stirling engine. 

The multi-objective optimization has been applied by different 
engineering problems such as skyline computation and vehicle 
routing issues [10-12]. Evolutionary algorithms (EA) were applied 
in an attempt to stochastically solve problems of this generic class 
during the 18th century [13], Inquiring a set of solutions, of which 
each satisfies the objectives at an acceptable level without being 
overcomed by any other solution, is a reasonable solving method 
to a multi-objective quandary [14], Multi-objective optimization 
of various thermodynamic and energy systems have been of high 
interest attention of researchers these days [15-21], 

In this study, the multi-objective optimization is conducted 
with four decision variables including the temperature of the heat 
source, the engine stroke, the mean effective pressure, and the 
engine frequency in order to increase the efficiency and power 
and reduce the pressure drop. Here, both second-order and third- 
order optimizations are applied so the parameters were considered 
according to the following parameter package order: 


- Efficiency and pressure drop. 

- Power and pressure drop. 

- Power, efficiency and pressure drop. 

2. Stirling system 

Stirling engine is a type of hot air engine which can produce 
mechanical work using heat transfer between the hot and cold 
sources. In the Stirling engine working cycle, the external heat is 
transferred to the working fluid at the maximum temperature of 
the cycle and the heat is rejected from the working fluid at the 
lowest temperature of the cycle. 

When the working fluid passes through the regenerator during 
the movement from the hot side to the cold side of the engine, 
some heat is absorbed by the regenerator. This heat is stored in 
the regenerator, and rejected to the working fluid with minimal 
thermal dissipation while the working fluid passes from the cold 
space to the hot space. While displacer moves reciprocating, it 
forces the working fluid to a reciprocating movement by passing 
through the regenerator between the hot side and the cold side 
[ 22 ], 

Various types of Stirling engines include: Gamma type Stirling 
engine, Beta type Stirling engine and Alpha type Stirling engine 
[24-26]. All of them have the same thermodynamic cycle, but 
the mechanical mechanisms of them are different. The Beta Stirling 
engine is the oldest structure of Stirling engine (Fig. 1 ). These types 
of engines use the combination of power piston and displacer [27], 
The construction of this engine is in such a way that both pistons 
are located linearly in a cylinder. The displacer piston moves the 
working fluid among the heater, the regenerator, and the cooler 
between the hot and cold spaces [28], The studied engine for this 
research was the GPU-3 whose characterizations are demonstrated 
in Table 1. 

3. Optimization 

Due to the necessity to optimize several objectives simulta¬ 
neously and access to several conflicting criteria, commonly the 
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Power piston 


Cooler 
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Displacer piston 


Cartei- 

Rod 

Crank 


Fig. 1. Beta type Stirling engine [23], 


Table 1 

GPU-3 engine specification [23], 


Specifications 

Values 

Type 

Beta 

Heat source T H 

977 K 

Heat sink T c 

288 K 

Mean pressure 

4.13 MPa 

Working fluid 

Helium 

Working frequency 

41.7 Hz 

Regenerator length 

0.0226 m 

Wire diameter 

40 rim 

Regenerator diameter 

0.0226 m 

Total mass of gas 

0.001135 kg 

Porosity 

0.697 

Displacer diameter 

0.069 

Number of regenerators per cylinder 

8 

Duct number of gauzes of the matrix 

308 


Non-dominated crowding distance 

sorting sorting P t+1 



Fig. 2. The structure of NSGA-11 multi-objective optimization model [32], 


One of the ways that has shown its effectiveness in solving 
many multi-objective problems, and has been used for 
multi-objective optimization algorithm, is the NSGA-II optimiza¬ 
tion algorithm [30], Fig. 2 shows the performance of NSGA-II 
multi-objective optimization model. 


3.1. Decision making 


For choosing the optimum point of the engine, a decision should 
be taken and there are different ways to make decisions. In this 
work, three methods including TOPSIS, LINMAP, and Fuzzy have 
been used which are explained briefly as followings. In the LINMAP 
method, the optimal point is a point on the Pareto Front that has 
the shortest distance to the ideal point. The ideal point is a point 
that has best situation. In the TOPSIS method, the optimal point 
is a point which has the lowest value of C L . 


where d,+ and d r are the distance of each point from the positive 
ideal point and negative ideal point, respectively. In the Fuzzy 
method, the Fuzzy Bellman-Zadeh is used which is explained in 
further details in the literature [33-35]. 


combination of goals are used to find out the solution of multi¬ 
objective problems, while these problems must be solved as 
multi-objective. Unlike the single-objective optimization which 
searches unit optimal solution; in the multi-objective problems, 
due to a conflict between the goals, there is no single optimal solu¬ 
tion. Thus, the method should search several optimal solutions that 
are called the Pareto set [29], 

There are several techniques for solving the multi-objective 
optimization problems, among which the Epsilon-Constraint 
Method, the Weighted Sum Method, and the Evolutionary 
Algorithm can be counted. In this study, the multi-objective evolu¬ 
tionary algorithm has been applied due to the ability to solve com¬ 
plex problems and provide the optimal exchange curves between 
the objectives. These models can easily solve the problems that 
do not adherence certain continuity, the decision spaces are not 
integrated, and their objective functions have random parameters. 
This model was proposed by Deb to solve the problems of the clas¬ 
sic genetic algorithm [30], The main problem of the previous 
multi-objective optimization includes a bulk of calculations per 
iteration that leads to increasing the run-time of modeling and 
do not maintain the appropriate number of optimal values during 
the implementation of the model [31], 


3.1.1. Bellman-Zadeh decision-making method 

As the intersection of all fuzzy criteria and constraints, a final 
decision is defined by the Bellman and Zadeh model and this model 
is represented by its membership function. In this method, a 
matrix of membership function is assigned, where the fuzzy mem¬ 
bership function is placed in the column for each objective. The 
rows of the membership matrix reflect the magnitudes of the 
membership functions for each solution of the Pareto frontier. 
Hence, the number of rows is equal to the number of the solutions 
and the number of columns in the membership function matrix is 
equal to the number of the objectives on the Pareto frontier. In the 
literature, the details of the method definition are described [33- 
35], For the proposed solution, the fuzzy membership function of 
each solution is set at the minimum value of the membership func¬ 
tions of all objectives. Therefore, a vector of membership is 
obtained in which its elements are the minimum membership 
function of the objectives at each solution. As the final solution, 
the minimum distance with the ideal point is selected for the best 
solution. In other word, in the Bellman-Zadeh fuzzy decision mak¬ 
ing method, a maximum of minimums (minimum membership 
function of objectives) is selected as the final optimal solution. 
For more clarifications, the details are reported in other works 
[33-35], 






























































































524 


S. Toghyani et al./Energy Conversion and Management 87 (2014) 521-529 


3.2. Objective functions, decision variables and constraints 


5. Theoretical model 


A mathematical model of the optimization process has three 
general parts (the objective functions, the decision variables and 
the constraints problem) that are described here in. Three objective 
functions are used for this study including the pressure drop inside 
the engine, the power and the efficiency of the Stirling engine. In 
this paper, four decision variables have been considered as 
followings: 

/: Engine’s frequency. 

P: Mean effective pressure. 

d: Stroke. 

T h : Temperature of heat source. 


In order to do conduct the optimization steps, the Nlog analyzer 
was applied while the code was modified with some complemen¬ 
tary codes for the linkage between the analyzer and the optimiza¬ 
tion codes. Nlog is a third-order Stirling cycle analyzer which 
calculates the generated power and rejected heat by Stirling engine. 

The code converts the whole gas transferring channels and 
pipes of the engine, to some divided control volumes and deter¬ 
mines the gas dynamic and thermodynamic parameters for each 
control volume by solving the continuity, momentum and energy 
equations simultaneously. The standard forms of continuity, 
momentum and energy equations are presented by Eqs. (6)-(8), 
respectively [37]. 


The objective functions with respect to the following 


constraints have been solved: 


20 s; / < 50 Hz 

(2) 

2 < P < 5 MPa 

(3) 

0.02 < d ^ 0.08 m 

(4) 

800 s: T„ 1100K 

(5) 

4. Nlog analyzer 



f s n,dA, + iy^<,v) 

l pv '" M ‘ + s{{ pdv ) 


= 0 


= 0 


( 6 ) 

(7) 


d_ 

dt 


^ j pVfdV^J + J PndAf + 


I, 


( v f p) v f n dA f + t fk A Wk 


= 0 


Nlog is a third order code in which, the engine is considered as a 
complex unit and the engine is divided into several parts. In order 
to achieve the heat transfer equations, the equations of mass, 
momentum and energy of the engine are written to describe the 
heat transfer processes for each part. In this method, the basic 
equations, in forms of partial differential equations, are written 
versus time and place. Then, some simplifications such as consid¬ 
ering one-dimensional conditions for each part and steady flow 
are carried out [36], Pressure drop and mechanical friction in the 
compression and expansion spaces are considered as losses. The 
following variables are entered into the Nlog code: 

- The Stirling engine type. 

- The geometric properties of all gas transferring channels, pipes, 
and chambers. 

- The wetted wall surface areas and roughness factor for each. 

- Regenerator’s type, geometry and porosity. 

- The properties of the non-ideal working gas (including com¬ 
pressibility factor, fluid thermal conductivity and viscosity 
relations). 

- The engine initial charge pressure and local temperatures. 

- The heat exchangers wall temperature gradient (constant along 
time period). 

- Stiffness of springs if connected to any moving parts, neglecting 
spring mass. 

- Geometries of the linkages between the engine moving parts (If 
exist). 

- The external time varying force/torque applied to the engine. 
The Nlog output consists of: 

- The engine output power and efficiency. 

- Pressure against total volume. 

- Heat transfer rate against time. 

- The working gas velocity in pipes and channels. 

- Rigid bodies’ velocity versus time. 

- The temperature of working gas at any point of the engine ver¬ 
sus time. 



The free flow area for the working gas transformation is assumed to 
be constant for each control volume, thus Eq. (6) is converted to Eq. 
(9). The area is also assumed to be constant during time, so Eq. (9) 
may be converted to Eq. (10) with this assumption. 


j s PVfndAf + ^- t (p k Af k x k ) = 0 
j s PPfndA f +A fk ^(p k x k ) = 0 


(9) 

( 10 ) 


The density and velocity of the working gas are assumed to be 
uniform over the cross sectional area between the adjacent control 
volumes and inside each control volume. Thus Eq. (11) becomes the 
final simplified continuity equation used in the code. 

Vf s M/i + A fk gt (p k x k ) = 0 (11) 

s=S k 

Extending the above assumptions of working gas velocity, density 
and cross sectional area, for working gas pressure in Eq. (7) will 
result in Eq. (12). The following is the simplified form of the 
momentum equation. 

§- t MtWfk) + + t/Ak = o (12) 

s=S k s—S k 


The energy equation is also simplified upon assuming the gas 
specific heat to be equal for every place of each control volume 
and is expressed as followings: 


( p A** ( Uk + -f )) + E ( p '*“* + Ps-j -+ p 0 v fMf< 


-AfP^-^O 
lk k dt dt 


( 13 ) 
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The total heat transferred into each control volume can be calcu 
lated using the following relation: 

dQc k _ Nu k K k A Wlc 
dt 


“CLc k 1 vU /< 1V k n w k / x T x 
~HT= d k (Tw k ~T k ) 


(14) 


The variation of wall temperature for each control volume is 
calculated from: 
dT Wl , 1 fdQ_ Ck dQe^ 
dt 


1 w k 


C w m„ k 


dt 


= 0 


(15) 


The gas specific internal energy is assumed to be only a function of 
gas temperature according to Eq. (16) for helium [38]: 


u k = CJ k ,(C„ = 3.11 kj/kg K) 


(16) 


The gas thermal conductivity is assumed to be a function of gas 
temperature and pressure according to Eq. (17) for helium [38], 


K k = 2.68 x 10 


~ 3 x (1 + 1,123P fc l0 s 


..7l(l-2xlO~ 9 P k ) 


(17) 


The thermodynamic ideal gas behavior is considered and is defined 
as followings [37], 

P k = p k rT k , (r = 2.08 kj/kg I<) (18) 

Two assumptions are used by Nlog code depending on the wall-type 
of each Stirling engine chambers. The assumptions for an isentropic 
regenerative type wall are shown in Eqs. (19.a) and (19.b) 
represents the assumption for the isothermal type. 

dt 
dT w 


= 0 


4 ^ = 0 

dt 


(19.a) 

(19-b) 


The formulations of Reynolds number (Re) and kinetic Reynolds 
number (Re w ) are described in Eqs. (20)-(23). 

Pk v fJ>K 


Re k = - 


Re Wk — 


P k D 2 h k 8 

Pk 


/.i k = 1.475 x 10 


80 


D/i,, = - 


44, 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


A schematic of the engine linkage kinematic parameters that are 
modeled by the code is shown in Fig. 3. 


- (e - r) 2 


b 2 =d{l- r) - e 2 


b e = \JV - (e + r) 

V e = V c ,e +A d (b 0 — b 2 — r sin 8) 
V e = V dc + 2 A p (b, - b e ) 


(24) 

(25) 

(26) 

(27) 

(28) 



Thus, V e to V c could be expressed as the functions of crank angle 
according to Eqs. (27) and (28) of which the parameters b u - b 2 , 
and b 0 are obtained from Eqs. (24)-(26), respectively. For each 
control volume which is characterized by k index, Eqs. (11)—(18), 
contain 12 unknown variables including x, t, Nu, v s , p, P, u, CL, CL, 
K, T w and T. The values of r and Nu are expressed as the functions 
of Reynolds number Re and kinetic Reynolds number Re w according 
to the experimental correlations which have been presented by 
Kays and London [39], de Monte et al. [40,41] and Zheng [42] for 
different ranges of Re, Re w different types and dimensions of cross 
sectional area, the wetted surface roughness can be obtained for 
the engine heat exchangers. 

Eqs. (11)—(19) totally form a system of nine equations which are 
solved by Nlog code to determine the nine remained unknown 
variables at each time step for each control volume. Finally, the 
code running results including the engine output power and the 
rejected heat are calculated using Eqs. (29) and (30). 


P = 


E ArA 

k:Va riable 
volumes 


dxk 

dt 


Qr = E 0* 

fe:Cooler 

volumes 


(29) 

(30) 


6. Results and discussion 

For the purpose of initial investigation on the final model for 
determining the objective and criteria, a preliminary study and 
modeling was carried out, which is firstly described as below. 



Fig. 4. Variation of the output power of the Stirling engine for different frequency at 
different strokes. 



Fig. 5. Variation of the efficiency of the Stirling engine for different frequency at 
different strokes. 
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Fig. 6. Variation of the pressure drop of the Stirling engine for different frequency 
at different strokes. 
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Fig. 9. Variation of the pressure drop of the Stirling engine for different frequency 
at different temperatures. 



Frequency (Hz) 


60 


Fig. 7. Variation of the output power of the Stirling engine for different frequency at 
different 515 temperatures. 



Fig. 8. Variation of the efficiency of the Stirling engine for different frequency at 
different temperatures. 


6.1. The effects of frequency, stroke and heat source temperature 

According to Fig. 3, the output power is increased up to a certain 
point by increasing frequency, after that, the power is decreased 
with increasing frequency. Also in certain frequency, with increas¬ 
ing the stroke, the output power is decreased. As it is demonstrated 
in Fig. 4, by increasing the frequency, the efficiency is decreased. 
Because, by increasing frequency, the working fluid has less time 
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Fig. 10. (a,b) Pareto optimal Frontier in the objectives space (pressure drop- 
efficiency). 

for heat exchanging between the hot source and the cold source. 
As it is shown in Fig. 5, by increasing frequency, the differential 
pressure drop is increased. It is happened due to creating pressure 
drop for increasing the engine fluctuations. All these results are 
reported for a constant temperature (see Fig 6). 

The effect of heat source temperature on the output power in 
different frequency is demonstrated in Fig. 7. As can be seen in 
Fig. 7, by changing the frequency, the output power has a maxi¬ 
mum power, and until a certain point the power is increased then 
decreased. Also in a certain frequency, by increasing heat source 
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Table 2 

The coefficients and goodness of fit for pressure drop-efficiency curve. 


Coefficient 
a = 3.629e+13 


b = 44.65 

c = 0.09088 

Goodness of fit 

SSE: 0.003777 

R-square: 0.9353 

Adjusted R-square: 0.9304 

RMSE: 0.01205 



& 

-O 
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• DPI vs. Power 
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Pressure Drop = axexp(bxPower) + 
cxexp(dx Power) 

R2 = 0.98 








4 


& 




(h) 


Power (kW) 



Fig. 12. (a and b) Pareto optimal Frontier in the objectives space (pressure drop- 
power-efficiency). 


Fig. 11. (a.b) Pareto optimal Frontier in the objectives space (pressure drop-power). 

temperature, the output power is increased. According to Fig. 8, the 
efficiency is decreased by increasing frequency, because the 
differential pressure drop is increased by increasing frequency 
(Fig. 9) and the working fluid has less time for heat exchanging 
between the hot source and the cold source. The results of this sec¬ 
tion are presented at constant stroke, 0.03, which is taken from the 
experimental work. 

6.2. The results of optimization 

For reaching the aim of the multi-objective optimization, the 
pressure drop (bar) of the Stirling engine must be minimized and 
the efficiency and the output power (kW) of the Stirling engine 


must be maximized. The optimization was carried out using a class 
of multi-objective evolutionary algorithm called NSGA-II. The 
optimization process was performed using the third order thermo¬ 
dynamic analysis with the constraints which were expressed in 
(2)-(5). Fig. 10a illustrates the Pareto Frontier in the space for pres¬ 
sure drop and efficiency. Three final solutions were selected by the 
Fuzzy, L1NMAP and TOPSIS decision-making methods which are 
indicated in this figure. In Fig. 10a, it is shown that the optimal 
values of the differential pressure are varied from 0.08 bar to 
0.28 bar and the optimal value of the efficiency is between 0.38 
and 0.48. The equation of the Pareto optimal Frontier is obtained 
and indicated in Fig. 10b, which can be applied to gain the system 
pressure drop, as followings: 

dP = a x if + c (31) 


Table 3 

The coefficients and goodness of fit for pressure drop-power curve. 


Coefficient 
a = -53.46 

b = 0.2064 

c = 53.49 

d = 0.2068 


Goodness of fit 

SSE: 0.02932 

R-square: 0.98 

Adjusted R-square: 0.9791 


RMSE: 0.02108 
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Table 4 

The coefficients and goodness of fit for pressure drop-power-efficiency curve. 

Coefficient 
pOO = -0.09471 
p02 = -14.24 

Goodness of fit 

SSE: 0.01437 R-square: 0.992 Adjusted R-square: 0.9908 RMSE: 0.01574 


pl0 =-0.4763 p01= 3.673 p20 = 0.02009 pll =2.182 

p30 = 0.001209 p21 =0.03438 pl2 =-2.848 p03 = 16.14 


Table 5 

Comparison of the final optimal solutions of efficiency-pressure drop (multi objective-multi variable) with experimental results. 


Methods 

Decision variables 



Objective functions 


MK) 

d (m) 

P(bar) 

/(Hz) 

Efficiency (%) 

Pressure drop (bar) 

Initial theoretical status 

977 

0.03 

41.3 

47.1 

36.17 

0.475 

TOPSIS 

90187 

0.0635 

21.25 

21.17 

43 

0.097 

LINMAP 

90,187 

0.0635 

21.25 

21.17 

43 

0.097 

Fuzzy 

902.96 

0.0631 

21.69 

21.34 

45 

0.1007 

Experiment [43] 

977 

0.03 

41.3 

47.1 

35 

- 


Table 6 

Comparison of the final optimal solutions of power-pressure drop (multi objective-multi variable) with experimental results. 


Methods 

Decision variables 



Objective functions 


MK) 

d (m) 

P (bar) 

/(Hz) 

Power (kW) 

Pressure drop (bar) 

Initial theoretical status 

977 

0.03 

41.3 

47.1 

4.17 

0.475 

TOPSIS 

1085.87 

0.0273 

50 

32 

4.316 

0.297 

LINMAP 

1083.87 

0.0243 

50 

31.73 

4.46 

0.325 

Fuzzy 

1022.89 

0.0253 

48.1 

3626 

4.7 

0.3988 

Experiment [43] 

977 

0.03 

41.3 

47.1 

3.9 

- 


Table 7 

Comparison of the final optimal solutions of power-efficiency-pressure drop (multi objective-multi variable) with experiment results. 


Methods 

Decision variables 



Objective functions 



Th (K) 

d (m) 

P (bar) 

/(Hz) 

Power (kW) 

Efficiency 

Pressure drop (bar) 

Initial theoretical status 

977 

0.03 

41.3 

47.1 

4.17 

36.17 

0.475 

TOPSIS 

1000.59 

0.029 

48.36 

34.93 

4.26 

42 

0.358 

LINMAP 

1000.59 

0.029 

48.36 

34.93 

4.26 

42 

0.358 

Fuzzy 

1022.62 

0.045 

47.1 

41.66 

4.51 

40.7 

0.447 

Experiment [43] 

977 

0.03 

41.3 

47.1 

3.9 

35 

- 


Table 8 

Comparison of the final optimal solutions of this paper (multi objective-multi variable) with corresponding results. 


Type of model 

/(Hz) 

P (bar) 

d (m) 

TH (K) 

Power (Kw) 

Efficiency (%) 

Initial theoretical model 

47.1 

41.3 

0.03 

977 

4.17 

36.17 

LINMAP and TOP SIS 

34.93 

48.36 

0.029 

1000.59 

4.26 

42 

Fuzzy 

41.66 

47.1 

0.045 

1022.62 

4.51 

40.7 

Urieli and Berchowitz (adiabatic model) [44] 

47A 

41.3 

0.03 

977 

8.3 

62.5 

Timoumi dynamic model without losses [43] 

47A 

41.3 

0.03 

977 

7.1 

54.96 

Timoumi dynamic model with loss dissipation (Ml) [45] 

47A 

41.3 

0.03 

977 

7.4 

53.1 

Urieli and Berchowitz quasi-steady flow (pressure drop included) [46] 

47A 

41.3 

0.03 

977 

6.7 

52.5 

Timoumi dynamic best model [46] 

47A 

41.3 

0.03 

977 

4.2 

38.49 

Experimental [43] 

47A 

41.3 

0.03 

977 

3.9 

35 

Error between initial theoretical model and experimental 

- 




6.9 

3.3 


The italic values of f- frequency; d- diameter; TH- hot side temperature. 

The coefficients (with 95% confidence bounds) and goodness of fit 
are summarized for the pressure drop and the efficiency in Table 2. 

The Pareto optimal Frontier for the output power and pressure 
drop is illustrated in Fig. 11a, that is a two objective function. It is 
clear that the optimal value of the differential pressure is varied 
from 0.1 bar to 0.68 bar, and the optimal value of the power is 
between 1 and 6.5 kW.The equation of the Pareto optimal Frontier 
is obtained and indicated in Fig. lib. The coefficients (with 95% 
confidence bounds) and goodness of fit are summarized for the 
pressure drop and the efficiency in Table 3. As it can be seen, if 


the objective functions would be power and pressure drop, the 
heat source temperatures for engine using LINMAP, TOPSIS, and 
FUZZY methods will be 1083.78, 1085.35, and 1022.894, respec¬ 
tively. Also, the strokes are obtained as 0.0273, 0.0243, 0.0253 
(m) using LINMAP, TOPSIS, and FUZZY methods, respectively. 

dP = a x exp(fa x p) + c x exp(d x p) (32) 

Pareto optimal Frontier for three objective functions, power, and 
efficiency of Stirling and system pressure drop are shown in 
Fig. 12a. The optimal results for the decision parameters and 
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objective functions using LINMAP, T0PS1S and Fuzzy decision¬ 
making methods are summarized in Table 4. 

The equation of the Pareto optimal Frontier is obtained and 
indicated in Fig. 12b, which can be applied to gain the system 
pressure losses, as followings: 

dp = p 00 + plO x p + pOl x i] + p20 x p A 2 + pi 1 x p x i) 

+ p02 x tj A 2 + p30 x pr A 3 + p21 x p A 2 x l] + pl2 
x power x p A 2 + p03 x p A 3 (33) 

where p is output power. 

The comparison of the optimal results obtained in this paper 
using three decision making methods with the experimental 
results of GPU-3 Stirling engine is demonstrated in Tables 5 and 

6. As it can be seen in Table 7, the output power and the efficiency 
of the engine are increased to 600 W and 5.7%, respectively, when 
using Fuzzy decision making. Also by the use of TOPSIS and LIN¬ 
MAP decision making, the output power and the efficiency of the 
engine are increased to 460 W and 7%, respectively. 

According to Table 8, the characteristics of this work are com¬ 
pared to the other works which have not considered the important 
losses. This reason caused that they reported higher efficiency. 

7. Conclusions 

In this study, the third order thermodynamic analysis has been 
considered to optimize the output power, the efficiency and the 
pressure drop of a Stirling engine. The efficiency-pressure drop, 
the power-pressure drop and the efficiency-power-pressure drop, 
of the Stirling engine were applied simultaneously for multi-objec¬ 
tive optimization and the temperature of heat source, the stroke, 
the engine pressure charge, and the engine frequency were consid¬ 
ered as the design parameters. It was observed that the Fuzzy deci¬ 
sion-making method yields in a final solution with more desirable 
output power and efficiency. Also, in aspect of pressure drop, the 
TOPSIS and LINMAP with the same amount have a better situation 
than the Fuzzy model. 
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